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In these notes I explain how to describe one-dimensional quantum systems that are simulta- 
neously near to, but not exactly at, a critical point, and in a far-from-equilibrium steady state. 
This description uses a density matrix on scattering states (of the type of Hershfield's density 
matrix), or equivalently a Gibbs-like ensemble of scattering states. The steady state I am con- 
sidering is one where there is a steady flow of energy along the chain, coming from the steady 
draining / filling of two far-away reservoirs put at different temperatures. The context I am 
using is that of massive relativistic quantum field theory, which is the framework for describing 
the region near quantum critical points in any universality class with translation invariance and 
with dynamical exponent z equal to 1. I show, in this completely general setup, that a particular 
steady-state density matrix occurs naturally from the physically motivated Keldysh formulation 
of nonequilibrium steady states. In this formulation the steady state occurs as a result of a large- 
time evolution from an initial state where two halves of the system are separately thermalized. 
I also show how this suggests a particular dependence (a "factorization" ) of the average current 
on the left and right temperatures. The idea of this density matrix was proposed already in a 
recent publication with my collaborator Denis Bernard, where we studied it in the context of 
conformal field theory. 
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1 Introduction 



These notes will explain how to describe quantum systems that are simultaneously near to (but 
not exactly at) a quantum critical point, and in certain steady states away from equilibrium. This 
situation is particularly intriguing because it mixes in an essential fashion both the rich collective 
quantum effects and nonequilibrium physics. Also, steady states, being time independent, are 
as near as possible to the usual states of equilibrium physics while still being out of equilibrium, 
and we can provide general descriptions in parallel to parts of equilibrium statistical physics. 

We have in mind quantum systems composed of a macroscopic number of microscopic con- 
stituents with local (few-neighbors) interactions on some lattice. In certain situations, for in- 
stance when a parameter (a magnetic field, an interaction strength, etc.) is adjusted appro- 
priately, the system becomes critical. Near to such a quantum critical points, quantum effects 
dominate the low-energy and large-distance physics. Indeed, on the one hand, the ground state 
encodes quantum correlations between local observables that exist much beyond few lattice sites. 
This is nontrivial, because with local interactions, one expects that constituents separated by 
many sites do not "feel" each other very much, which is indeed what happens away from critical 
points. On the other hand, at least intuitively, the low-lying excitations can be understood as 
large-scale collective behaviors of the underlying constituents that emerge thanks to quantum 
fluctuations; as if large groups started to act together forming new quantum objects. Critical 
points arguably give some of the most interesting physics: they are where many nontrivial phys- 
ical effects arise (one of the most well-known being the Kondo effect); and their large-distance / 
low-energy physics is actually universal, because most of the intricacies of the local interactions 
have been lost in the averaging out over many sites that arises in large-distance correlations or 
collective objects. 

The physical theory that describes what is observed "near" critical points is quantum field 
theory (QFT) (more precisely, when not exactly at the critical point, massive QFT). By "near" 
we mean in the limit where observation lengths, lengths related to the system's geometry and 
correlation lengths are all large in lattice spacings, these three kind of lengths being taken to 
infinity in fixed proportion^]. This is the scaling limit. QFT is the physical theory that will 
provide predictions for what happens in the scaling limit (for instance, how correlations decay 
to zero, giving both the exponents and universal amplitudes). Of course, in reality one never 
reaches exactly the scaling limit. But when criticality is approached in macroscopic systems, 
then every large-distance and low-energy observations are effectively predicted by QFT. 

Note that if the correlation lengths are much larger than everything else, then we are exactly 
"at" the critical point. In these cases, QFT possesses scale invariance, and often with this, in- 

1 The concept of observation length is kept slightly nebulous here, but as examples this could be the distance 
between two points in a correlation function; or, translating into the energy-momentum language by the usual 
Fourier-transform inversion, the temperature. 
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variance under the full global conformal group. In 1+1 dimensions (the scaling limit of quantum 
chains), this very often comes with local conformal invariance, and we have a conformal field 
theory (CFT). We will not describe these situations here, and refer the reader to [1] for the CFT 
counterpart of the nonequilibrium steady states we will study here. 

The equilibrium physics of critical and near-critical points has been very well studied over 
the years, and there is by now a wide range of techniques to obtain exact and approximate 
results. Not as much can be said of the physics away from equilibrium. One reason is, of course, 
the large variety of out-of-equilibrium situations (instead of saying what the situation is, we say 
what it is not...). But also, even in the relatively simple cases where a steady state exists, there 
lacks a framework as powerful as and as widely applicable as that of Gibbs ensembles. Hence, 
one often must revert to other formulations, usually rather adapted to the type of physical 
model at hand. Here, we will start, in the massive QFT context, with an "explicit", real-time 
construction of the steady state, where we implement theoretically the idealized experimental 
set-up that should produce the steady stat^E This will lead us to rather general statements, 
paralleling Gibbs ensembles. 



2 Basic notions 



Consider a quantum chain of length L with local (nearest-neighbor, say) interactions. We want 
to describe a state where energy flows along the chain due to a difference of temperatures at its 
extremities. There are various ways to do this: one could use Lindblad's theory of open quantum 
systems where the thermal baths have been traced out and the time evolution is non-unitary [2] , 
or, in a somewhat related fashion, one could explicitly connect Caldeira-Leggett baths [3] with 
different temperatures on the extremities of the chain (as is done, for instance, in [U [5]). One 
also could consider an infinite chain, instead of a finite one, and try to describe the flow using 
asymptotic states in a Gibbs-like ensemble, through what is usually referred to in the condensed 
matter community as Hershfield's density matrix (Hershfield considered this for general charge 
transport [6]; Ruelle obtained rigorous general results in the context of C*-algebras with similar 
underlying ideas for heat transport [7], see also [8]); this of course requires a clear principle to fix 
this density matrix. Finally, also with boundaries sent to infinity, one could explicitly construct 
the baths with different temperatures by using parts of the system itself as baths: the infinite 
chain is initially cut into two semi-infinite parts, each thermalized at different temperatures, 
which are then connected and evolved unitarily for a long time. This latter method is in a sense 
the most "physical", as the steady state is explicitly constructed in real time (it is sometimes 
referred to as the "partitioning approach" of Caroli et al. [9], and sometimes as the Keldysh 
formulation, used for instance in the context of mesoscopic systems jlOl , and it is related to 
quantum quenches). Note that in this method, thanks to the unitary evolution, no part of the 
system is kept at a fixed temperature during time evolution; the information of the temperature 
is only in the initial mixed state. When combined with ideas of perturbation theory and path 
integrals, the Keldysh formulation leads to the powerful Keldysh perturbation theory (in fact, 
Keldysh perturbation theory can be used as well with the method of the Cladeira-Leggett baths). 

In the present notes, we will start with the Keldysh formulation in the context of general, 
nonperturbative massive QFT, and we will show, using non-rigorous but well motivated QFT 

2 Similar real-time constructions, where unitary evolution is performed from a state that is not a Hamiltonian 
eigenstate, are usually referred to, in general, as "quenches." 
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arguments, that it leads to a Hershfield formulation, giving the precise Gibbs-like ensemble of 
asymptotic states describing the steady state where there is an energy flow. This was already 
proposed in pQ, where it was studied in the context of CFT (there, the heat current and its full 
counting statistics were found) . See also [HI Q21 [131 EH CG2 EE E3 CEB1 CE] for similar constructions 
concerning charge or heat transport in a variety of other situations and at various levels of rigor. 

To be more specific, we consider the Hamiltonians H l and W of the two semi-infinite, 
separate halves (left and right) of a quantum chain, and the initial density matrix 

where the halves are at different temperatures T^ r = /3f r - We then evolve this density matrix 
with the full Hamiltonian H = H l + H r + AH where AH represents the energy term for the 
single link connecting both halves: 

pit) = e- iHt p e iHt . (2.2) 

The steady state average of observables is obtained upon evaluating lim t _ s . 00 Tr (p(t)0) /Tr (p(t)). 
More precisely, if v is a typical velocity of quantum propagation in the chain (which we will take 
below to be 1), then the limit is that vt must be much greater than the observation length 
(which could be taken as the length of the support of the observable O), and much smaller 
than the length of the system itself (which we have implicitly taken as infinity). This is so 
that both halves of the system are big enough to play the role of thermal baths, absorbing and 
emitting energy as necessary, so that a nonequilibrium steady state may arise. See [12|, [T4"l [T7] 
for discussions about such steady-state limits. 



Remark 2.1 One may worry about the physical meaning, in the Hershfield and Keldysh formu- 
lations, of sending the extremities of the chain to infinity. Indeed, according to Fourier's law, 
the heat current (which is, in many situations, the same as the energy current) is proportional 
to the gradient of the "local temperature" : j = — cVT with heat conductivity c. But if the ex- 
tremities are infinitely far from each other, and no part of the system is actually kept at a fixed 
temperature during the evolution, surely the local temperature gradient, after infinite time, will 
be zero, hence there should be no heat current, unless the heat conductivity c is infinite. Yet one 
of course find that the heat current is nonzero (although we will not evaluate any current in the 
present notes). We must realize that in low dimensions, Fourier's law is known to be broken: 
heat conduction is anomalous 11 9v \2(A [21\/ . and gets non-local contributions. In particular, it was 
shown in [28] that a nonzero Drude weight plays an important role. An underlying principle in 
the present context may be that of quantum coherence. Returning to Fourier's law for the sake 
of the argument, it should be understood in a coarse-grained fashion: one must take the limit of 
T(x + Ax) /Ax where Ax is much smaller than the distance at which T changes significantly, 
but much greater than the distance at which things start to mix and scatter (so that, in partic- 
ular, the concept of a local temperature is clear). But quantum coherence precludes scattering 
and mixing; so Ax must be much greater than the quantum coherence length. Near to quantum 
criticality, the correlation length is very large (in lattice spacings), hence one may indeed expect 
the quantum transport behavior to be very different from that predicted by Fourier's law. 

Remark 2.2 We will consider the case where there is a quantum critical point, and take the 
scaling limit, described by QFT. See J7?| [22^ for a discussion of the interplay between the scaling 
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limit and the steady-state limit. In particular, the scaling limit describes the low-energy behaviors, 
and one may worry that some low-energy states of H l + H r are high energy states of H and are 
not captured by a QFT description of the steady-state limit (i.e. lie far outside the low-energy, 
universal range). The basic explanation \22 f is that if the renormalization of the observable 
with respect to the H l + H r ground state is also finite with respect to the H ground state, then 
the correct steady-state average of this observable is captured by QFT. 

Remark 2.3 Instead of connecting a single link, we could also think of connecting a whole finite- 
size (in lattice spacing) subpart of the chain (i.e. AH corresponds to a finite number of sites). 
Near to critical points, this will lead to the same low-energy results, because the correlation length 
is infinite in lattice spacings. In fact, we also expect that with AH corresponding to a number 
of site being a finite proportion of the correlation length, we still find the same steady state. 

3 Energy-flow steady states in massive 1+1-dimensional quan- 
tum field theory 

We restrict our attention to translation-invariant quantum chains in the thermodynamic limit 
and assume that there is a quantum critical point with dynamical exponent z = 1. Averages 
near to, or at, quantum critical points have asymptotic behaviors described by QFT. With 
z = 1, the near-critical (not exactly critical) behavior of the quantum chain is described by 
massive relativistic QFT, which is the context that we consider. See for instance the standard 
reference [23j for a description of quantum critical points. In essence, if m is the gap between 
the lowest energy state of the quantum chain and the next state, and J is some microscopic 
energy scale (like the intersite coupling), then the scaling limit leading to massive QFT is that 
where T/, T r <, ~ m <C J. 

3.1 Asymptotic states and field configurations 

One way of thinking about QFT is as a theory for quantized field. Although this is exactly what 
the name relates to, it is not entirely accurate because in particular of the need for renormal- 
ization; yet it does provide a good intuitive picture, and gives rise to precise results when used 
correctly. In this picture, QFT is a quantum theory whose Hilbert space C has for basis the 
set {|</>)} parametrized by all field configurations <fi on real space, with overlap (4>\4>') = 6[<f>,(jJ] 
the delta-functional supported on the field configurations being equall. In the 1+1-dimensional 
case, these are field configurations <fi : R — > S on the line, where the field takes values in some 
space S. There is the associated x-dependent operator 4>{x) (which, following tradition, we refer 
to as a field) defined by <j)(x)\4>) = 4>(x)\<p), and there is an energy density h(x), formed out of 
<j)(x) and its derivatives, using which the Hamiltonian is 

/oo 
dxh(x). (3.1) 
-oo 

We can think intuitively of a field configuration as being the conglomeration of the individual 
basis states that the microscopic components of the underlying quantum model are in; but we 

3 It is not clear that this is actually a Hilbert space: restriction to the set of field configurations may be 
necessary, beside other subtleties; but for our intuitive description, this is sufficient. 
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must realize that on any finite length on IR, there are infinitely many such components, and that 
in the effective description of QFT, the space S is not necessarily the set of such microscopic- 
component basis states. 

A somewhat surprising consequence of the principles of QFT (locality and Poincare invari- 
ance) is that the Hamiltonian eigenfunctions and eigenvalues are completely characterized by 
few parameters: the particle spectrum. Indeed, given a set of quantum numbers AT and associ- 
ated particle masses A4 = {m a : a £ J\f}, a set of orthonormal basis states in the Hilbert space 
is given by 

{\0i,---,O k ) ai ,..., ak :0 l >...>6 k , aj eM, k = 0,1,2,...} (3.2) 
on which the Hamiltonian has eigenvalues 

k 

m a . cosh 9j . (3-3) 

3=1 

Here we are using the rapidity 6, in terms of which the relativistic (on-shell) energy and momen- 
tum are given by E = m cosh 6 and p = m sinh 6. When k = 0, we denote the state by |vac). The 
sets Af and A4 are characteristics of the QFT model, although they do not completely define it. 

Remark 3.1 Hence, in massive relativistic QFT, the problem of diagonalizing the Hamiltonian 
is almost solved: one only needs to determine few parameters! Determining the masses can 
be hard of course, and even when this is done, the QFT model is not yet solved, neither even 
fully defined. One more piece of information is necessary: that having to do with locality. One 
needs to say what the stress-energy tensor is, as a (tensorial) operator. This not only provides a 
definition of space via the momentum operator, but also of what operators are local: one can write 
the quantum locality condition of vanishing commutators with the energy density at space-like 
distances. This fully defines the QFT model, and to solve it (according to a certain definition), 
one needs to find all local operators and calculate their correlation functions. This in principle 
provides the measurable physical information encoded in the QFT model. 

The basis (I3.2p is a basis of asymptotic states. Naturally, asymptotic states are particular 
linear combinations of field configurations, formally 

■ ■ ,h)a u ...,cx k = ^^,...,^(01, ...,6k \<t>) \4>)- (3-4) 

Intuitively, they are linear combinations such that if we were to evolve the state, say, back in 
time for a very long time, we would find large wave packets widely separated and propagating 
at rapidities 9%,..., Of.- The wave packets would be associated to the various particle types 
a\, . . . so they would be wave packets of the various "fundamental fields" related to the 
various particle types. These various fundamental fields can be components of 4>(x) (remember 
that (j) takes values in S, which may be W 1 for a n-component field for instance), or some 
composite fields (like 4>{x) 2 , appropriately renormalized) , depending on the model. 

Of course, the paragraph above cannot be entirely true, because asymptotic states diagonalize 
the Hamiltonian; hence they evolve trivially (only getting a phase). A more correct statement 
is that one must construct linear combinations with the properties of the above paragraph, and 
then take the limit of infinite time, simultaneously infinitely-large and infinitely-separated wave 
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packets. In order to represent scattering states, this limit must be taken in such a way that 
the "free" relativistic evolution of wave packets would produce space-time trajectories (for the 
centers of the wave packets) that are found, at time 0, in a finite region of space. Where exactly 
they are found (so-called impact factors) is up to us, and just affect our choice of basis elements. 

Remark 3.2 What is described here are in-states. There is a similar construction where the 
condition is obtained by evolving forward in time; the resulting asymptotic states are out-states. 
These then provide two bases for the Hilbert space that diagonalize the energy operator; they 
have the same description and lead to the same energies, but are different bases. The overlaps 
between in-states and out-states forms the scattering matrix of the model. It is expected that 
giving the scattering matrix is equivalent to specifying the stress-energy tensor: it provides the 
locality information of the QFT model. This is rather explicit in the context of integrable QFT, 
with the so-called form factor program, but not obvious out of integrability. 

Remark 3.3 Not all possible linear combinations ^^c^c/)) of field configurations occur in the 
basis ( Iff. 2\) (or infinite linear combinations of basis elements); in particular the states \4>) them- 
selves do not occur. But we expect that by taking the completion, involving all (possibly infinite, 
and continuously so) linear combinations that converge, one recovers an appropriate set of field 
configurations \<j>) and their (again possibly infinite, and continuously so) linear combinations 
Yl(f> c 4>\ ( i ) ) with well-behaved coefficients c<f,. In fact, we expect to obtain all square-integrable 
field configurations (j) - this is natural, as <j) is supposed to represent a wave function in the 
"second-quantization" viewpoint on QFT. This is similar to the situation in ordinary quantum 
mechanics, where the wave functions for energy eigenstates may have certain continuity and dif- 
ferentiability properties, and appropriate infinite linear combinations need to be taken in order 
to obtain all possible square-integrable functions, including for instance the discontinuous ones. 
In the present case, however, this is much more complicated, because the basis \3. 2\) has contin- 
uous components (it is not a countable basis); hence the full Hilbert space structure (including 
completeness) is much more subtle. We will avoid these tricky details here, and only discuss 
basic properties in due course. 

It is extremely important to realize that asymptotic states are just particular linear com- 
binations of field configurations. The term "asymptotic" refers here to the way these linear 
combinations are chosen: by an asymptotic condition on a hypothetical time evolution. It does 
not refer to the states "living only in the far past" (or the far future). It turns out that this con- 
struction produces linear combinations of field configurations that diagonalize the Hamiltonian. 

We can verify the latter sentence and make the intuitive construction above much more 
precise as follows. Recall first that the vacuum in quantum field theory is itself a nontrivial 
linear combination of field configurations: 

|vac) = c (vac |0) \<j>). (3.5) 

One standard way of defining this linear combination is via the path integral (appropriately 
renormalized). Assume that the Hamiltonian is associated to the action functional S[ip], for 
ip : M. 2 — > S, where tp(x,t) is the trajectory in time t of the field configuration (f>. From this it is 
of course possible to write down the action Sji[ip] on any open region R of R 2 . One defines the 
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coefficients c(vac\(ft) by 



c(vac|0) = / D^e^M (3.6) 

tp(x,t)^0 : igl, i->-(l-iO+)oo 

where L = £ I 2 : t < 0} and is an integration measure over all configurations of 

(p : L — > S. That is, one performs Feynman's path integral on field trajectories over half of 
space-time, the half with negative times, with the boundary condition in the limit of zero time 
t — > 0~ given by the field configuration (ft, and with infinite-past asymptotic condition (under 
Feynman's imaginary-time prescription) equal to (here we assume that (ft = is a minimum 
of the potential representing the interaction). 

The wave packets are constructed on top of this vacuum. Let us denote by (ft a (x) a funda- 
mental field associated to the particle a. This means that (ft a (x) has the quantum numbers of 
the particle a, and that the Fourier transform of the two-point correlation function 

(vac|T (cft a (x, t)(ft a (0, 0)^ |vac) (3.7) 

(where T is the time-ordering operation, bringing the earliest operator to the right, and (ft a (x, t) = 
e lHt (ft a (x)e~ lHt ) has a pole, as function of the square two-momentum E 2 — p 2 , at rr? a and no 
other singularities at lower values. We assume for simplicity of the discussion that (ft a is a real 
bosonic field, transforming under the Poincare group in the spin-0 representation (generalizing 
to other representations and statistics and to complex fields is straightforward, if somewhat 
tedious). The statement above implies that the correlation functions of (ft a (x,t) with any other 
local fields satisfy the Klein-Gordon equation asymptotically in time (other equations occur for 
other representations): 

m2 -^-§2)( T {^( x ^---)) =0 { e ~ m,y/TI ^) as t 2 -x 2 ^oo. (3.8) 

Here m! is the lowest mass, greater than m, created by (ft a , and • • • represents other fields at 
fixed positions. 

Define the following operators: 

f°° ( d - d \ 

A a {9) = i hrn^ J dx yf e (x,t)—(f) a (x,t) - -^fe{x,t)(ft a (x,t) J (3.9) 



with 

(x - coth(6)ty 



fe(x,t) = exp 



im cosh((9) t — im sinh((9) x 



L- 



(3.10) 



as well as their Hermitian conjugate Aa(9). Here the limit on t and L is such that t ^> L, so 
that the wave packets are well separated (without going into more details). Using (|3.8p and 
taking into account the limits taken in (|3.9p . one can show that these are eigenoperators of the 
Hamiltonian, in particular: 

[H,Ai(6)\ = m a cosh 8 Al(9). (3.11) 

The vacuum having zero energy, this means that we can construct eigenstates of the Hamiltonian 
by multiple actions of Aa(9): 



n, 



h) au ...,a k =Al 1 (0 1 )-..Ai k (9 k )\v*c). (3.12) 



S 



These are the asymptotic (in-)states. 

One may wonder why we didn't use the operators A a {9) instead in order to create different 
Hamiltonian eigenstates out of the vacuum. This is because 

A a (6)\v&c) = (3.13) 

and 

[A a (e),Al,(9')} = ^6^6(9 - 9'). (3.14) 

The former is a (slightly subtle) consequence of (|3.5j) , (|3.6p and (13. 9|) , (I3.10|) . The latter is more 
involved, and a consequence of locality and of the pole structure of the Fourier transforms of 
the two-point functions (vac|T (^(f> a (x, t)(f> a /(0, 0)^ |vac). Relations f)3. 13|) and 13.14]) mean that 
the vacuum (|3.5p can alternatively be defined as the vacuum of the Fock space of the canonical 
algebra generated by A a (9) and Aa(9); this is a definition which does not rely on the existence 
of an action functional. 

Remark 3.4 Note that the canonical commutation relations j3.14\ ) do not imply that the model 
is free. They only relate to the fact that the asymptotic states have a particularly simple descrip- 
tion, which can be obtained as the Fock space of a canonical algebra. But for interacting models, 
the local fields are not simple Fourier transforms of the operators A a (9), Aa(9), and taking the 
limit t — > oo instead of t — >■ — oo in 113. 9\) gives rise to operators creating the out-states, which 
are in general different form the instates. 

Out of the states (|3.12p . we may now construct the Hilbert space of the model by completion. 
Without going into any detail of this rather subtle mathematical process (see Remark I3.3f) , we 
will denote by H the resulting Hilbert space; this is expected to be a proper subspace of the 
space of field configurations C. 

3.2 Asymptotic states in the presence of a cut-impurity 

Impurities and boundaries on the quantum chain give rise to QFT models with inhomogeneities 
running, in space-time, along the time direction (hence there is a natural, preferred time direc- 
tion). These can be treated in similar ways as above. The case of interest below is that where 
there is a cut at one point along the chain (that is, we have two separate half- infinite chains). In 
this case, the field configurations 4> are on M° := K\ {0} (the basis \<f>) forms the space C°), and 
the Hamiltonian H° is (slightly) different from H, in its expression in terms of fields </>, having 
a discontinuity at x = 0: 

H° = f dxh(x). (3.15) 

Obviously, we can write H° = H l + H r as a sum of a term on (— oo, 0) and a term on (0, oo), 
respectively. These terms commute with each other, [H l ,H r ] = 0, hence can be diagonalized 
simultaneously 

The vacuum can be described as 

|vac)° = c°(vac|0) \<j>) (3.16) 
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where 

c°(^) = J ^ _^ x):xeu0 Vuxp&M (3.17) 

with L° = {(x,t) £ R 2 : x ^ 0,t < 0}. Here we choose the free boundary conditions as x 0~ 
and x — > + , hence we do not constrain the limit value of the field. This is clearly a factorized 
vacuum, |vac)° = (vac)' (g> |vac) r . 

Then, asymptotic states are constructed in the same way as for the model on the whole line. 
We define A° a by 

f'°° ( d - 8 \ 

A° a (e) = ilm lo / dx(fe(x,t)^l(x 3 t)- m f e (x,t)^(x 3 t)) (3.18) 

with (f>a(x,t) = e lHOt (j) a {x)e~ tHOt . Then we simply act on the vacuum with the operators 

1*1. - , hf ai ,..., ak = «) f (0i) • • • (A° a J(0 k )\v & c)°. (3.19) 

The l-r factorization also occurs for asymptotic states. Note that for positive rapidities, the 
wave packets (|3.10p are supported on (— oo, 0), and for negative rapidities, on (0, oo) (apart 
from exponentially decaying tails which vanish in the limit t — > —oo). Since asymptotic state 
operators commute for rapidities of different signs, we can then write 

\Ou---,0 k )° ai ,..., ak = II (O^lvac/ ® J] (^.)%)|vacr. (3.20) 
i:6i>0 j:0j<0 

The factors are the separate eigenstates of H l and H r . 

Again, it is important to recall that the asymptotic states (I3.19P are just particular linear 
combinations of field configurations on R°. In particular, the field configurations involved, 
seen as functions on K, are generically discontinuous at x = (that is, the left- and right- 
limits generically don't agree). Despite the similar formulation, the asymptotic states (|3. 19[) 
are very different from the asymptotic states (]3.12p . when seen as linear combinations of field 
configurations (in particular, there is reflection on the point x = in the case of the theory on 
E°). 



Remark 3.5 A model on MP is equivalent to two independent copies of the model on a half-line. 
In each copy, the set of asymptotic states is twice smaller than that for the model on the line, 
and one can describe them by choosing a sign of the rapidities. This is clear from the asymptotic 
state construction above: for instance, only states where particles have positive momenta have 
wave packets supported inside (— oo,0), as noted. Putting two copies together, we get a set 
of asymptotic states isomorphic to that of the model on the line. There are different ways of 
describing the resulting tensor product states. But if one copy is interpreted as being on (— oo, 0) 
and the other on (0, oo) (so that we have a model on R°, as we did above), then we have both 
signs of rapidities involved, and the asymptotic state construction that we have gives the same 
description h3. 2\) both for the set of asymptotic states on R, and that on R°. Here we neglect 
states with particles of zero rapidities, because they contribute only infinitesimally to averages of 
local operators. 
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Remark 3.6 Impurities and boundaries give rise to more intricate scattering phenomena, which 
include in general reflection and transmission phenomena. This is not seen in the construction 
of asymptotic instates, or of asymptotic out-states; it is rather seen in the structure of the 
scattering matrix, formed by the overlaps between in- and out-states. The scattering matrix will 
have additional components corresponding to reflections and transmissions. In the case of the 
cut-impurity, of course, only reflection phenomena will occur at x = 0, without transmission. 
In general, one could also have additional bound states to the impurity; but this is not expected 
to occur with a cut impurity. 

Out of the states (|3.19p . we may now construct the Hilbert space 7i° = H 1 ®% r of the model 
by completion. We omit any discussion of the subtleties of how rl° compares to W (e.g. how it is 
embedded into H), but note that the two Hilbert spaces are closely related, since only one point 
on the line is missing in the construction of Prom the viewpoint of the asymptotic state 
basis, the only difference between H and H° is of course the states with zero-rapidity particles, 
which, as mentioned, are expected to contribute only infinitesimally to averages of local (or 
finite-extent) observables. 

Remark 3.7 It is not essential to single out the point x = as we have done above, since in 
QFT, the value of the field at one point in space does not affect physical results (unless infinite 
potentials exist at that point). From this viewpoint, we may expect H = H, and we could simply 
add the point at x = from a configuration on MP in some more or less arbitrary fashion (e.g.: 
the average of the right and left limits). However, singling out x = clarifies the discussion and 
avoids these technical details. We may also think of the point x = to be representing the scaling 
limit of some finite, contiguous-sites subsystem of the underlying quantum chain that separates 
the half-infinite left and right parts. Our results apply to this general situation. In fact, we may 
replace the point x = by any finite interval [a, b], a, b € M. On the left / right of the interval, 
we may use the description of % l ® Ti r as above, and we would have an additional space £ to 
describe the configurations on the interval itself. Our main results would still hold in this case. 



3.3 Steady-state density matrix for energy flows 
3.3.1 Statement of the result 

According to the steady-state construction in Section [21 and taking the scaling limit, we need to 
evaluate the average of observables H -evolved to time t, in the mixed state with density matrix 
po given by (12. ID . The left half of the system, with energy H l , is at temperature j37 , and the 
right half, with energy H r , is at temperature fi~ x . The steady state is obtained in the limit where 
t — > oo. For our observable, we will take any operator O that is finitely- supported. Formally, 
this is an operator that commutes with the energy density h(x) whenever \x\ is large enough; 
informally, it is composed of local fields (in products, integrals) at points that all lie on some 
finite closed interval (the support). This guarantees that far enough in space, the observable 
has no effect. 

We will show that 

ljm T rw ,[ P0 ^ Oe -"'] _ !*„ [p„„01 
t^co Tv n o[p \ Tr^[p sta tJ 
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for every finitely-supported operator O. Here, po is diagonal on the asymptotic-state basis (|3.19p 
of rl° with eigenvalues 



(3.22) 



where E l and E r are the eigenvalues of H l and H r , respectively; and p s tat is diagonal on the 
asymptotic-state basis (|3.12p of % with eigenvalues 



-p t E+-p r E-. 



(3.23) 
(3.24) 



where 

E + = E l = m ai cosh6>,, E~ = E r = m aj coshOj. 

i:0i>O j:8j<0 

Note that E + and E~ have the same form as E l and E r , respectively, when expressed in terms of 
the rapidities of asymptotic states. However, we emphasize that the asymptotic states associated 
to the values E^ are very different from those associated to the values E l ' r . 

Using the density operators n a {9) = Aa(9)A a (9), which measure the density of asymptotic 
particles of type a at rapidity 9, we can write explicitly p sta t as 



Pstat = exp 



-Pi 



/•oo rO 

/ d9j2 n *( e )-Pr d#J> a (0) 
JO J-oo 



(3.25) 



Hence we have a description of the non-equilibrium steady state via a density matrix; a "quasi- 
equilibrium" description. 



Remark 3.8 Note that E + +E~ is the eigenvalue of H on the state $3.12\) . but that {3iE + +(3 r E~ 
is not, in general, the eigenvalue of a local charge (the integration of a local density). This is in 
contrast to the situation on TiP : not only E l + E r is the eigenvalue of H° , but also j3iE + + j3 r E~ 
is the eigenvalue of fiiH 1 + (3 r H r , a local charge. Nonlocality of the charge involved in the 
nonequilibriun steady-state density matrix is a quite general phenomenon, observed also in the 
cases of charge currents in impurity models 241- 



Remark 3.9 Although we discuss density matrices and traces over them, there are a great 
many subtleties involved into a proper definition of such objects. One way to define them is 
to regularize to a finite-length quantum chain in the same universality class, and then take the 
limit after evaluating the trace. This is correct but in general rather messy, and the difficulty 
in taking the scaling limit makes it hard to extract the universal structures associated to QFT 
from the particularities of the quantum chain chosen. Another way is to consider the operator 
algebra of QFT independently of the states (C* -algebra), and to see the trace as a particular 
linear functional on this algebra. This is the algebraic- QFT way, is mush more elegant, and can 
be directly applied to the infinite-length limit, which is useful given that this limit must be taken 
first in the construction of the steady state. Unfortunately, it is hard to do anything else than 
the free theory (and aspects of conformal field theory) in this completely rigorous picture. Taking 
ideas from both perspectives may be the most fruitful way. 

The derivation of (|3.2ip is based on the statement that the operator e - lHt e lHOt ; j n the large-t 
limit, is a map T~L° — > % which takes the state (I3.19p . times any state on C, to the state (I3.12P : 

hm e-^V^^!,...,^)^,...^ = A) Q1 ,...,a fc . (3.26) 
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This statement is true in the context of evaluating matrix elements of finitely-supported oper- 
ators (that is, in the "weak" sense); in particular, there is no problem with a possible loss of 
invertibility after the limit process (see Remark l3.7p . This statement shows (|3.2ip . because, H° 
commuting with po, we have 



Ti- 



ft 



[po- 



AHt 



Oe 



-iHtl 



Tr 



H° 



iHt„-iH°tQ e -iHt e iH°t 



Tr H o [po 



Tr w [po 



(3.27) 



and we just have to take the trace by summing over the states (|3.19p . and use the statement on 
every term (see Remark 13.91 for the subtleties in principle involved in taking the trace) . 



Remark 3.10 The result derived below can be seen as naturally related to the so-called G ell- 
Mann and Low theorem of QFT 125]. This theorem essentially implies that if the limit on the 
left-hand side of \3.26\) exist, then it must give the right-hand side. There are, however, many 
subtleties in the existence of this limit, and of the limit on the trace $3.21\) itself; in particular, 
the theorem itself involves an adiabatic process, where as here we have a quench. The arguments 
we present are based on the particular setup of the Keldysh formulation, with the particular 
order of limits involved, different from what is used in standard QFT constructions. Hence we 
believe it is worth going through the derivation without the explicit use of the Gell-mann and 
Low theorem. 



Remark 3.11 In integrable systems, one expects that time evolutions from states that are not 
Hamiltonian eigenstates (quenches) give rise, after a large time, to Generalized Gibbs Ensem- 
bles 126\/ . These ensembles are described by density matrices of the form Pqucnch — ^ J ^ ^ 
where Qj are the local conserved charges (including the Hamitonian) and where the generalized 
(inverse) temperatures (3j can be determined by requiring that the averages of the Qj before the 
start of the time evolution be reproduced by p qU ench • In the cases of integrable QFT, it is tempting 
to view our density matrix /o s tat ^ n this context. However, our result holds beyond integrability. 
Further, contrary to the usual quench situation, the density matrix p s t a t does not reproduce the 
conserved averages of the charges Qj. This is because /3 sta t reproduces averages of local fields 
and describes the steady-state region only, while conserved charges Qj are integration of local 
densities over the whole length of space including the far-away baths. The subtlety lies in the 
order of the limits that we must take: infinite length first, then infinite time. 



3.3.2 Existence of steady state 

As a first step, we will show that the steady state exists, and that, in particular, we can replace 
the limit t — > oo by the Wick rotated limit, 

lim , (3.28) 

t=— IT, T—>00 

for the operator e -" lHt e lH ° t that's on the right of O in (13.271) (and t = ir for the one on the left). 
In order to show this, we write 



e -iHt e iH°t = Vexp 



-t 

(3.29) 



i / dnSH^K) 
. Jo 
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where SH = H - H° = dx /i(0) and SH(k) = e m ° K 5H e~ iH ° K . Then, we have 



Tr 



p e -iH°t e iH 



Q e -iHt e iH°t 



Tr 



po V exp 



; j°_ t d,K 8H(k) 



O Texp 



; J * dn 5H{t 



Tr [po' 



-iH t f ,iHt ( ,-iHt f ,iH t'\ 



Tr 

"P exp 



Po exp 
o 



V exp 



6H(k) 



O Pexp 



J * dK,5H(n 
-t 

dn 5H(k 



Then, we use the fact that for local fields in QFT, we have the form 
<Oi(0)O 2 (k)>S 



connected _ K #g-#TM g im K ag i | 



OO 



connected 

(3.30) 

o 



(3.31) 



and similar equations for time separation of groups of fields. We use this with the local fields 
being O and various products of 5H(k) at various times k. This implies that the integrals in 
(|3.30p all converge as t — > oo, hence existence of the steady state. Further, in the path-ordered 
exponential on the right we can replace t *— > —it (and then in the integral k = ik) and we get 
convergent integrals still (opposite shift on the left), with the same value; hence we can use 



3.3.3 Vacuum mapping 

Since the asymptotic states f)3. 19|) are constructed from the vacuum |vac)°, we start by describing 
the effect of e ~ lHt e lH 1 Q n |vac)°. We employ a standard argument from elementary QFT, used 
in textbooks in order to obtain, for instance, the perturbative series from the operator formalism. 

Let us write |vac)° as a linear combination of eigenstates of H: 

|vac)° = B vac |vac) + B *\ s ) = -Svaclvac) + J de"^ B a {6)\6) a + ... 

excited states s 

Then 



e- im |vac) 



| vac) + B s e 

excited states s 



-iE.t 



We have shown the large-t limit can be taken by doing (I3.28|) . Then all contributions from the 
excited states vanish, and we get 



lim e -iHt e iHHu m „\0 



t—>co 



I vac) 



B v 



vac) 



(3.32) 



Hence, in the large-time limit, the vacuum |vac)° collapses to |vac), as long as the overlap i? vac 
is nonzero (in which case it will normalize away by taking the ratio of traces in (|3.2ip ). 

In order to justify _B vac ^ 0, we use arguments similar to those of standard QFT textbooks, 
except for some changes due to the construction of the steady state requiring us to take the 
length of the system infinite first; the arguments then use in an essential way the locality of 
5H = H - H°. 

In standard QFT, in the context of perturbation theory, one compares a "true" vacuum 
| vac) with the vacuum of a free theory |vac)°. The differences between the Hamiltonian is a 
term integrated over the whole volume - here the length of the system I - hence the overlap 
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(vac|vac)° ~ e decreases exponentially with I. Hence one usually assumes that the volume is 
finite, and take the limit of infinite volume afterwards; this gives rise to the usual perturbation 
theory. If we had a similar situation here, this would be very bad, since for the construction 
of the steady state, we need to take the limit I — > oo before t — > oo. Happily, the difference 
between H° and H is a local term (in fact, it is an infinitesimal term, but this is not the main 
point). This implies that the overlap B V . AC = (vac|vac)° is stable in the limit £ — > oo. Indeed, in 
the language of path integrals, we can represent it as 

(vac | vac) 

^(vacjvac) °(vac|vac)° / 'j V&<pe iS fV R o xR <pe iS 

where H = {(x,t) : x S R, t > 0}. Let us divide space-time into cells of extent much larger 
than the inverse of the smallest mass, mT 1 (but still finite). Then the large-volume asymptotics 
of the path integral is obtained by factorizing it into such cells. Using time-reversal symmetry, 
we see that every cell, except three, cancel out in the ratio; in particular, for every cell, in 
the numerator, containing a space-time region (t < 0, x « 0) cut into two pieces, there are two 
equal cells, in the denominator under the square root sign, containing a similar space-time region 
it < 0, x ~ and t > 0, x ~ 0). Those not cancelled are those containing the space-time point 
(0,0), because they are different in the three overlaps involved. Since only three cells remain, 
the result has a finite infinite- volume limit. Here we do not discuss UV divergencies: these are 
present, but one can UV regularize the theory and take the full scaling limit only at the end of 
the calculation of the nonequilibrium average (see Remarks 12.21 and I3.9P . 



/ Z>L0uHye tS 



3.3.4 Excited states mapping 

Finally, in order to map the full excited states as in (|3.26[) . we only need to show that (see (|3,9p 
and (pTIHl) ) 

lim e- im e iHOt A° a {6)e- iHOt e iHt = A a {6). (3.33) 

t— >oo 

The main idea is to equate the "asymptotic-state" time t = t asymv used, in (|3.9p and (]3.18p . to 
construct the asymptotic states, with the negative —t steady of the "steady-state" time t = t steady, 
used for instance in (I3.33p . to construct the steady state. This should be fine by the fact that 
the observable O is of finite extent. Then (|3.33p follows simply from 

e- iHt e iHOt 4>° a (x, -t)e- iHH e iHt = 4> a (x, -t). (3.34) 

In order to justify equating t aS ymp = —tsteady, consider the trace (|3.27p . We may use, 
instead of the asymptotic states, a set of states created as in (13T9L but with A° a replaced by 
an expression like (|3.18p where the t limit has not been taken yet (hence neither the L limit), 
fixing t = t asymp . This set is not a basis of states, but it is likely that such states, at least 
for tasymp negative enough, span (a dense subset of) H . Hence the trace can be written as a 
sum over such almost-asymptotic states, with appropriate normalization coefficients. Denoting 
\s)t asyrnp the almost-asymptotic state corresponding to the asymptotic state |s)°, and C®{t asymp ) 
the associated coefficients (satisfying lim iasym;) _ > _ 00 C®(t asymp ) = 1), we have formally 

Tr w o[p •] «l> s °C s °(ia TO ) ta3 J p (s\ ■ \s)° tasymp (3.35) 
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where iv® is the weight (|3.22p . Denoting \s)t asymp the similar almost-asymptotic states corre- 
sponding to (|3.9p . and C s (t aS ymp), we have similarly 

Tr w [p stat • ] » w s C s (t asymp ) tasymp (s\ • \s)t asymp (3-36) 

s 

where w s is the weight (|3.23p . Now choosing t s t ea dy = —tasymp, Relation (|3.34p and the equality 
(|3.24p for the weights imply that we obtain, for the trace (13.271) . 

J2 W sC° s (-t s teady) -t ateady (s\0\s)-t 3teady • (3.37) 
s 

For general O, not necessarily of finite extent, we couldn't expect the limit t s t ea dy — )• oo of this 
expression to give the trace on the right-hand side of (|3.2ip , because it is the different coefficient 
C 9 (t steady) that would be required. However, for observables of finite extent, we expect 

, lim ^W^I^Wnp = (*\0\s). (3.38) 



Since the coefficient C®(—t s teady) itself tends to 1, the large steady-state time can be taken, and 
we obtain (13.2ip . 

Remark 3.12 Note that it is crucial to take the limit of infinite system length first, as otherwise 
we clearly cannot replace the asymptotic- state time by the (negative of the) steady-state time 
(with finite lengths, the large time limit won't give asymptotic states). 



3.4 A consequence on the average current 

The statement (I3.2ip with (I3.23P does not provide an immediate way of evaluating the average 
current or its fluctuations in general (in CFT it does, see the exact results in [T]). Yet, an 
intuitive argument leads to the nice observation of left- and right- "factorizaton" of the average 
current. Clearly, the average energy current is 

j = Tr^ [/) s tat p(0)] 
Tr w [pstat] 

where p(x) is the momentum density. In cases where there is a sense in the notion of quasi- 
particles and quasi-momenta even in finite volume (for instance, in integrable models via the 



Bethe ansatz, see for instance [27]), then there is a natural finite- length L extrapolation p^ t 
of 

Pstat- There is of course a finite- L notion of p, and by translation invariance, we have J 
lim^oo J L with 

-!^n L [pLtP] 



at 



Jl = L~ 



II LrstatJ 

where P = f_ L ^ 2 dx p(x) is the total momentum. Clearly, for every state \v) with n quasi- 
particles of momenta pj, j = 1, . . . , n, we have 
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This means that we can write P = P+ + P_ where P± are operators that measure the momenta 
of positive- and negative-momentum particles: 



right-moving particles, /) s tat = pstat,+pstat,-, similar to each other but depending on /3i and j3 r 
respectively, we obtain 



The above derivation holds even if the quasi-particle / quasi-momentum description is not exact 
at L finite, as long as it becomes "more and more exact" as L — > oo, which should be the case 
for any near-critical or critical system. 

The "factorization" relation (|3.39p . as well as an equivalent factorization for the moment 
generating function, was established in the context of CFT in pQ; further, interestingly, (13.391) 
was observed numerically in [28] in spin chain models, hence beyond QFT0. As was remarked 
in [28] . (|3.39p immediately implies a simple relation between the linear conductance and the 
nonequilibrium current: J = J(T h T T ) = ff r dT G(T) where G{T) = dJ{T' ,T)/dT\ T > =T . 

4 Concluding remarks 

We have shown how the Keldysh formulation leads, by QFT arguments, to a Gibbs-like ensemble 
formulation for nonequilibrium steady states of energy flow in quantum systems near to critical 
points. 

We note that all the characteristics of the constructions were essential in order to obtain the 
proof: that the system's length be infinite, otherwise we cannot identify asymptotic-state and 
steady-state times, and that the quench be local (the term 5H = H — H° is local), or at least 
of finite extent, otherwise B vac vanishes. The existence proof of the steady state presented also 
uses both the infinite length and the locality (or finite extend) of 5H. 

Trying to generalize this construction to higher dimensions lead to interesting problems. 
If we have one additional dimension of finite length, then everything goes through essentially 
unchanged (except for subtleties in defining asymptotic states). However, as can be expected 
physically, an additional dimension of infinite length will cause problems. There, one would 
expect, immediately after connection of the two halves, to have an infinite total flow of energy, 
while in the steady state itself, to have a finite total flow of energy hence a vanishing energy 
current density. This slightly "anomalous" situation is translated into the vanishing of the 
overlap B vac : in the cell-cancellation argument presented, there remains infinitely many cells, 
leading to a vanishing proportional to e~ £ where £ is the extent of the additional dimension. 

4 To my knowledge, Relation (|3.39|l was first presented in my talk about the work [T] with Denis Bernard at 
the Workshop on Non-equilibrium Physics and Asymmetric Exclusion Processes, ICMS, Edinburgh, 8 December 

2011. There it was proposed in the context of general 1+1-dimensional QFT. It was also presented in my talk 
at the Galileo Galilei Institute workshop "New quantum states of matter in and out of equilibrium" on 23 April 

2012, and appeared in the first version arXiv : 1202 . 0239vl of the preprint of Q]. 




J = f(f3l)-f(Pr). 



(3.39) 
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As we mentioned, the idea of certain Gibbs-like ensembles of asymptotic states representing 
a nonequilibrium steady state, and its link with the Keldysh formulation, has been used many 
times in the past. Hershfield's density matrix [6] was developed for charge transport, and was 
studied and used afterwards in the context of quantum impurities [321 ECU Ell EEH ESI E] . Charge 
transport through impurities has also been studied using the asymptotic states of integrable QFT 
[29 1 130 } I3T1 132] . and using Bethe-ansatz asymptotic states [33], without explicit reference to a 
density matrix (although it is in some sense there implicitly). Nonequilibrium Gibbs ensembles 
of asymptotic states in spin chains were discussed with mathematical precision in [3 \TT\ [8] 
for heat transport, and were also studied in [18] in free fermion and boson models. Energy 
transport, including exact results for the current and the so-called full counting statistics, was 
studied in the context of CFT in [I], where a nonequilibrium density matrix was obtained in 
terms of the Virasoro algebra. The present notes give the first general study of energy transport 
in the context of massive relativistic QFT, providing the density matrix for nonequilibrium heat 
flows (which was proposed in [1]) and a direct connection with the Keldysh formulation. 

Much still needs to be done, of course. The evaluation of the energy current and its fluctua- 
tions is foremost: can concepts of integr ability, combined with the present Gibbs-like ensemble 
description, be fruitfully used (I have work in progress with my collaborators on this subject)? 
Could one develop an efficient perturbation theory, including all the subtleties of renormaliza- 
tion, paralleling that for finite-temperature QFT? Deeper questions about the thermodynamic 
of nonequilibrium steady states may also be accessible from the present formulation. I hope that 
these notes will open the way to further developments in these directions. 



Acknowledgements 

I am very much indebted to my collaborator Denis Bernard for many (still ongoing) discus- 
sions about the subject of nonequilibrium steady states, which ultimately led to my understand- 
ing of the situation in massive QFT. I also thank the participants to the Galileo Galilei Institute 
workshop "New quantum states of matter in and out of equilibrium" for questions raised on the 
subject which eventually led to this work, and I acknowledge in particular discussions with J. 
Bhaseen and F. Essler. I am of course very grateful to the students and lecturers of this Les 
Houches summer school for their insightful questions and comments, and in particular to Natan 
Andrei for comments and encouragement to write these notes. 

References 

[1] D. Bernard and B. Doyon, Heat flow in non-equilibrium conformal field theory, J. Phys. 
A: Math. Theor. 45 (2012) 362001, preprint larXiv : 1202 . 02391 

[2] G. Lindblad, On the generators of quantum dynamical semigroups, Commun. Math. Phys. 
48 (1976) 119-130. 

[3] A.O. Caldeira and A.J. Leggett, Influence of Dissipation on Quantum Tunneling in Macro- 
scopic Systems, Phys. Rev. Lett. 46, 211 (1981). 

[4] K. Saito and A. Dhar, Fluctuation theorem in quantum heat conduction, Phys. Rev. Lett. 
99 (2007) 180601, preprint [arXiv : cond-mat/0703777l 



18 



[5] K. Saito and A. Dhar, Generating Function Formula of Heat Transfer in Harmonic Net- 
works, Phys. Rev. E 83, 041121 (2011), preprint ErXiv : 1012 . 06221 

[6] S. Hershfield, Reformulation of steady state nonequilibrium quantum statistical mechanics, 
Phys. Rev. Lett. 70, 2134 (1993). 

[7] D. Ruelle, Natural nonequilibrium states in quantum statistical mechanics, J. Stat. Phys. 
98, 57 (2000), preprint [arXiv : math-ph/9906005| 

[8] W.H. Aschbacher, V. Jaksic, Y. Pautrat and C.-A. Pillet, Topics in non-equilibrium sta- 
tistical mechanics, in Open Quantum Systems III, Springer, pp. 1-66, Aug. 2006, Lecture 
Notes in Mathematics. 

[9] C. Caroli, R. Combescot, P. Nozieres and D. Saint-James, Direct calculation of the tun- 
neling current, J. Phys. C 4, 916 (1971). 

[10] A. -P. Jauho, N. S. Wingreen and Y. Meir., Time-dependent transport in interacting 
and noninteracting resonant-tunneling systems, Phys. Rev. B 50, 5528 (1994), preprint 
larXiv : cond-mat/9404027l 

[11] W.H. Aschbacher and C.-A. Pillet, Non-equilibrium steady states of the XY chain, J. Stat. 
Phys. 112 (2003) 1153. 

[12] B. Doyon and N. Andrei, Universal aspects of non-equilibrium currents in a quantum dot, 
Phys. Rev. B 73 245326 (2006), preprint larXiv : cond-mat/0506235l 

[13] J.E. Han, Mapping of strongly correlated steady-state nonequilibrium system to an effec- 
tive equilibrium, Phys. Rev. B 75, 125122 (2007), preprint [arXiv : cond-mat/0604583| 

[14] B. Doyon, The density matrix for quantum impurities out of equilibrium, lec- 
ture notes for the Fifth Capri Spring School on Transport in Nanostructure, 
|http : //tf pi . physik . uni-f reiburg . de/Capri09. 

[15] P. Dutt, J. Koch, J. Han, K. Le Hur, Effective equilibrium theory of nonequilibrium 
quantum transport, Ann. Phys. 326 (2011) 2963-2999, preprint [arXiv : 1 101 . 15261 

[16] V. Moldoveanu, H.D. Cornean, C.-A. Pillet, Non-equilibrium steady-states for interacting 
open systems: exact results, Phys. Rev. B 84, 075464 (2011), preprint larXiv : 1104.5 399, 

[17] D. Bernard and B. Doyon, Full Counting Statistics in the Resonant-Level Model, J. Math. 
Phys. 53, 122302 (2012), preprint liTXiv : 1 105 . 16951 

[18] A. Dhar, K. Saito and P. Hanggi, Nonequilibrium density matrix description of steady 
state quantum transport, Phys. Rev. E 85, 011126 (2012), preprint ErXiv : 1106 . 32071 

[19] F. Bonetto, J.L. Lebowitz and L. Rey-Bellet, Fourier's Law: a Challenge for Theorists, in 
Mathematical Physics 2000, edited by A. Fokas et. al. (Imperial College Press, London, 
2000), p. 128, preprint [arXiv : math-ph/0002052| 

[20] S. Lepri, R. Livi and A. Politi, Thermal conduction in classical low-dimensional lattices, 
Phys. Rep. 377, 1 (2003), larXiv: cond-mat/01 12193] 

[21] A. Dhar, Heat Transport in low-dimensional systems, Adv. Phys. 57, 457 (2008), preprint 
larXiv : 0808 . 32"56l 



19 



[22] D. Bernard and B. Doyon, Nonequilibrium steady states in conformal field theory, in 
preparation. 

[23] S. Sachdev, Quantum Phase Transitions, Second Edition, Cambridge University Press, 
New York, 2011. 

[24] B. Doyon, New method for studying steady states in quantum impurity problems: 
the interacting resonant level model, Phys. Rev. Lett. 99 076806 (2007), preprint 
|cond-mat/0703249[ 

[25] M. Gell-Mann and F. Low, Bound states in quantum field theory, Phys. Rev. 84, 350 
(1951). 

[26] M. Rigol, V. Dunjko, V. Yurovsky and M. Olshanii, Relaxation in a completely inte- 
grable many-body quantum system: an ab initio study of the dynamics of the highly 
excited states of ID lattice hard-core bosons, Phys. Rev. Lett. 98, 050405 (2007), preprint 
larXiv : cond-mat/0604476l 

[27] L.D. Faddeev, How algebraic Bethe ansatz works for integrable models, in Les 
Houches 1995, Relativistic gravitation and gravitational radiation pp. 149-219, preprint 
|arXiv:hep-th/9605187] 1996. 

[28] C. Karrasch, R. Ilan, J.E. Moore, Nonequilibrium thermal transport and its relation to 
linear response, preprint arXi v: 1211 . 2236. 

[29] P. Fendley, A.W.W. Ludwig and H. Saleur, Exact non-equilibrium transport through point 
contacts in quantum wires and fractional quantum Hall devices, Phys. Rev. B52 (1995) 
8934, preprint [arXiv : cond-mat/9503172| 

[30] R. Konik, H. Saleur and A. Ludwig, Transport through Quantum Dots: Analytic Results 
from Integrability, Phys. Rev. Lett. 87, 236801 (2001), preprint |arXiv : cond-mat/0010270l 

[31] R. Konik, H. Saleur and A. Ludwig, Transport in Quantum Dots from the Integrability of 
the Anderson Model, Phys. Rev. B 66, 125304 (2002), preprint larXiv : cond-mat/0103044l 

[32] E. Boulat and H. Saleur, Exact low temperature results for transport properties 
of the interacting resonant level model, Phys. Rev. B 77, 033409 (2008), preprint 
larXiv : cond-mat/0703545l 

[33] P. Mehta and N. Andrei, Nonequilibrium Transport in Quantum Impurity Mod- 
els: The Bethe Ansatz for Open Systems, Phys. Rev. Lett. 96, 216802 (2006), 
preprint arXiv: cond-mat/0508026; P. Mehta, S.-P. Chao and N. Andrei, erratum, 
cond-mat/0703426 



20 



